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A detailed analysis of convergence properties of the Andreozzi-Lee-Suzuki iteration method, which 
is used for the calculation of low-momentum effective potentials Vi ow t is presented. After summa- 
rizing different modifications of the iteration method for one-flavor channel we introduce a simple 
model in order to study the generalization of the iteration method to multi-flavor channels. The 
failure of a straightforward generalization is discussed. The introduction of a channel- dependent 
cutoff cures the conceptual and technical problems. This novel method has already been applied 
successfully for realistic hyperon-nucleon interactions. 

PACS numbers: 13.75.Ev, 21.30.-x, ll.10.Hi 
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I. INTRODUCTION 

In Ref. pj Bogner et al. have developed a low- 
momentum NN interaction Vj ow k which is constructed 
by a combination of renormalization group (RG) ideas 
and effective field theory approaches. Starting from a re- 
alistic nucleon-nucleon interaction in momentum space 
for the two-nucleon system, the Vj ow k is obtained by 
integrating out the high-momentum components of the 
potential. In this way a cutoff momentum A is intro- 
duced to specify the boundary between the low- and 
high-momentum space. Since the physical low-energy 
quantities such as phase shifts or the deuteron binding 
energy must not depend on the cutoff, the full-on-shell 
T-matrix has to be preserved for the relative momenta 
k, k' < A. This yields an exact RG flow equation for a 
cutoff-dependent effective potential Vl ow k, which is based 
on a modified Lippmann-Schwingcr equation (LSE). 

Instead of solving the RG flow equation, which is a dif- 
ferential equation, with standard methods (e.g. Rungc- 
Kutta method) directly, an iteration method is used. 
This iteration method, pioneered by Andreozzi, Lee and 
Suzuki (ALS) @,H0, 

is based on a similarity transfor- 
mation and its solution corresponds to solving the flow 
equation ||.This transformation is also used for the ef- 
fective interaction in the no-core shell model (e.g. 0,0)- 
In Ref. H the Vi ow k approach was generalized to 
the hyperon-nucleon {YN) interaction with strangeness 
S = — 1. For hypcrons with strangeness S = — 1, two 
isospin states 7=1/2 and 1 = 3/2 are available. For 
/ = 1/2, several coupled hyperon-nucleon channels oc- 
cur that require new technical developments in the iter- 
ation method as already outlined in Ref. [8|- Especially, 
the AA" <-> SAT transition involves a mass difference of 
w 78McV and causes the failure of the known meth- 



ods. In the present study we discuss these novel technical 
modifications in detail and investigate the convergence of 
the ALS iteration methods. 

The paper is organized as follows: In the next Sec. we 
summarize the ALS iteration method for the one-flavor 
channel and discuss several modifications which improve 
its convergence. The generalization of the ALS method 
to multi-flavor channels is presented in the following Sec- 
tion. Instead of working with a realistic hyperon-nucleon 
interaction we introduce a simple solvable model in order 
to elucidate the essential modifications of the iteration 
method. We end with a summary and outlook. 



II. ONE-FLAVOR CHANNEL 

The so-called model space (or P-space) is introduced in 
order to truncate the infinitely many degrees of freedom 
of a Hilbert space to a finite, low-energy subspacc. The 
projection operator P onto the model space and its com- 
plement projection operator Q onto the excluded space 
arc defined by 



P 



J2\<t>i){<t>i 



Q = l-P 



(1) 



where d denotes the dimension of the model space. The 
wave functions \4>i) are eigenfunctions of the unperturbed 
Hamiltonian Hq, which is composed of the kinetic energy 
and an arbitrary one-body potential. 

Employing this separation of the Hilbert space, the 
original full Hamiltonian H can then be replaced by a 
transformed Hamiltonian H. that acts only on the low- 
energy model subspace. In particular, if j^j) is an eigen- 
state of the full Hamiltonian H then its projection onto 
the model space yi^f 1 ) = P\^i) satisfies the Schrodingcr 
equation 
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The crucial point is that the eigenvalues Ei of the 
transformed Hamiltonian are the same as those of the 
original Hamiltonian. 

On the other hand, the transformation of all model 
states V^) back into the exact states defines a non- 
singular wave operator via = f2 jV'f 7 )- Assuming 
that this wave operator has an inverse, the transformed 
Hamiltonian is related to the full Hamiltonian by a sim- 
ilarity transformation 



(3) 



There is no unique representation for the wave oper- 
ator il. Different choices for f2 may then give different 
results for the perturbativc expansion of the transformed 
interaction. Furthermore, even the same defining equa- 
tion for f2 can be solved by different iterative schemes. 
As an example one often employs the ansatz by Lee and 
Suzuki 0, Q for the wave operator 



(4) 



where u) is known as the correlation operator. The 
correlation operator generates the component of the wave 
function in the Q-space and can be expanded in terms of 
Feynmann-Goldstone diagrams [ij. One choice for the 
correlation operator is 



UJ 



QujP 



(5) 



and thus we have Q = 1 + uj. 

Applying the operator Q on Eq. ©, one immediately 
obtains QHP = 0. Together with Eq. Jsj, this then leads 
to the so-called decoupling equation 



luPHQuj + ujPHP - QHQuj - QHP = 



(G) 



which is to be solved in uj by iteration. For this non-linear 
matrix equation no general solution method is known, 
nor the number of existing solutions. 

In this work we will discuss and generalize the solution 
of the decoupling equation for arbitrary coupled chan- 
nels. The construction of the solution is based on the 
iteration method after Lee and Suzuki 0, E3 an d An- 
dreozzi Q, labeled as 'ALS iteration' in the following. 



A. The ALS iteration 

We start with an overview of the main ingredients of 
the standard ALS iteration. As already mentioned, this is 
not the only existing method which solves the decoupling 
equation. 

Defining the model space effective Hamiltonian, 

p (uj) = PHP = PHP + PHQuj , (7) 
and its complement Q-space Hamiltonian, 

q (uj) = QHQ = QHQ - ujPHQ , (8) 



an iterative solution lo so i = a of Eq. JBJ) can be obtained 
by 



cr = ^2 x n 



(9) 



where the s n 's satisfy successively 

-f 



Xq 



;QHP , 



QHQ 
xi = — — -x p (x ) 



l 



q (xo H h X n -l) 



X n ^ip (xq H 1" X n -l) 



Once a solution of these equations is known, the low- 
momentum effective Hamiltonian can be expressed as 



(10) 



Subtracting from 7Yi ow & the kinetic energy in the model 
space, PHqP, this yields the wanted low-momentum po- 
tential V\ ow fc. As a remark we note, that this iterative 
solution a is equivalent to the LS vertex renormalization 
solution (see Refs. 0,0). 

Since x n — > for large n, only a finite number of terms 
in the sum @ are needed numerically in order to reach 
a desired precision. Furthermore, this condition implies 
for each iteration step n the relation 



\\m\ < Wm\ i 



(11) 



where A p m is the largest P-space eigenvalue (absolute 
value) and X qm the smallest Q-space eigenvalue. Thus, as 
a necessary condition, this iteration method converges to 
the d eigenvalues of the Hamiltonian of smallest absolute 
value. 

A slight variation of the standard ALS iteration con- 
sists in the subtraction of a constant multiple of the iden- 
tity, mo/, from the full Hamiltonian H. In this case the 
iteration will converge to the shifted set of d eigenval- 
ues which are nearest to the chosen constant too. The 
choice of too is not arbitrary. One has to make sure 
that the eigenvalues closest to mo are really inside the 
P-space. For numerical applications, this variation can 
also be used to accelerate the iteration since the required 
shift of the eigenvalues is smaller if too is included. Fur- 
thermore, calculatin g th e V\ OV! k for the hyperon-nucleon 
Nijmegen potential |lH | (cf. Ref. @) the constant shift 
in the iteration is useful to avoid an unphysical bound 
state, which is inherent in the 1 Sq isospin 3/2 channel 
(cf. Ref. E3). 

The important observation is so far that each of the 
iteration schemes converges only to a set of d consecutive 
eigenvalues of the Hamiltonian, which arc supposed to be 
ordered in a nondecreasing order. 
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Andrcozzi described yet another generalization of the 
iteration procedure that enables one to shift the eigenval- 
ues individually by a set of different numbers instead of 
just one constant multiple of the identity (cf. Ref. [l3j K 
If one denotes by (A)i the i-th column of a matrix A, the 
decoupling equation for this case can be written in the 
form 

LoPHQ(Lo). t + lo (PHP - M) i - (QHQ - mi I q ) (lj) 1 
- (QHP) i = for i = 1, . . . , d (12) 

where the d-dimcnsional diagonal matrix M = 
diag (mi , . . . , m<j) contains the arbitrarily chosen num- 
bers rrii. Thus, for each column there is a separate de- 
coupling equation to be solved. Using the preceding it- 
eration procedure for the i-th equation one obtains the 
solution 



Mi 



1 



q (c„_i 



rriilq 



Xn-l^K-O-M); . (13) 



This modified ALS iteration scheme will converge to 
the eigenvalues closest to the chosen set of numbers rrii. 

This modification leads also to a stronger convergence 
condition for each iteration step n (cf. Eq. H11|0 



I A. 



(14) 



where X p ^ denotes the i-th eigenvalue in the model space 
and \ q j all eigenvalues in the Q-space. Thus, for a fixed 
eigenvalue i condition i|14|) there arc n — d inequalities for 
all eigenvalues in the Q-space. 

All iteration methods, presented so far, converge 
rapidly for channels with only one mass involved (one- 
flavor channel). For coupled channels with different 
masses (multi-flavor channels) new phenomena emerge 
and the convergence behavior (if any) of the iteration 
methods changes dramatically as already indicated in 
Ref. j8(. This becomes especially relevant e.g. for the 
hypcron-nuclcon interaction in the particle basis, where 
up to four particles with different masses are coupled. 



III. MULTI-FLAVOR CHANNELS 

As already mentioned, the standard and modified ALS 
iteration methods can fail for coupled channels interac- 
tions including different masses. In this Section we will 
assess the reasons for this failure in detail. Instead of con- 
sidering the realistic hyperon-nucleon interaction with 
different particles, we introduce a simple and solvable 
model which mimics all qualitative and relevant features 
of the realistic interaction. On the basis of this model 
the technical modifications of the ALS iteration method 
needed as a remedy for the coupled multi-flavor channels 
are discussed. 

In the chosen model the failure of the standard and 
modified ALS iteration methods for coupled multi-flavor 
channels and its cure are seen immediately We will re- 
strict ourselves to a coupled two-flavor channel in order 



to keep the discussion simple. The generalization to arbi- 
trary coupled multi-flavor channels with different masses 
is straightforward. 



A. A simple model 

The model is governed by the Hamiltonian H(k,k') = 
Ho(k) + V(k, k') which is in this case a 2 x 2 matrix in 
flavor-channel space. For the kinetic energy of the i-th 
channel we use 



1,2 



(15) 



with nil = 9 and m 2 = 13. All quantities in our 
model arc dimensionless. The reduced mass is defined 
by fit = mo?rii/(mo + rrii) with a fixed mass m = 8 
for each i-th channel. The mass parameters rrii are ad- 
justed in such a way that typical convergence problems of 
the iteration method, emerging in a realistic multi-flavor 
channel calculation, also appear in this model. 

As a bare Hermitian interaction we choose a Fourier- 
transformed Yukawa potential 



Vij{k,k') 



gjj l0 J Akk ' + b l 



*,j = l,2; (16) 



with the parameters an = 4.3, a\i = 2.0, 022 = 3.7 and 
6n = —3.0, 6 12 = —8.0, 6 2 2 = —4.0. The parameters 
dij reflect the depths and the bij the effective ranges of 
the corresponding potentials. They are chosen in such a 
way that the diagonal potentials Vu are larger than the 
off-diagonal ones. 

In the two-flavor channel, the block structure of the 
Hamiltonian looks like 



H 



_ fH u H 
H21 H 2 



12 



Hqa + Vu V12 
V21 Hq 9 + V22 



(17) 



with the relation V2i(fc, k') = V\2{k' \k). 

For each flavor-channel a P-space (Q-space) with two 
mesh points and a cutoff Ap = 5 (Aq = 10) is introduced. 
Here we use only two mesh points in order to simplify our 
discussion of convergence. In realistic calculations we use 
about 64 mesh points for each channel. 

Thus, each subblock Hij becomes now a 4 x 4 matrix 



Hi 



PHijP PHijQ 
QH l3 P QH 2j Q 



i,3 



1,2 



(18) 



In Fig. ^ the kinetic energy Eq. i|15|) as function of the 
momentum k for the i-th flavor channel is shown. Our 
choice of the P- and Q-spacc discretization yields four 
mesh points on each curve (squares for the first channel 
and bullets for the second channel). They are chosen in 
such a way that two mesh points of each flavor channel 
are close to the P-space cutoff Ap which is depicted by 
a vertical line in the figure. 



2 4 6 8 10 

k 

FIG. 1: The kinetic energies Ho,i as function of the mo- 
mentum k for the model. The line with squares denotes the 
first-flavor (i — 1) channel, the other one the second-flavor 
(i = 2) channel (see text for details). 

In the following we apply the standard and 
modified ALS iteration methods, described in the 
previous Section, to this model. The results are 
then compared with the exact solution of the 
model. It is solved by numerical diagonalization, 
which produces the following eight eigenvalues H = 
diag(9.55, 14.50, 15.05, 24.04, 13.64, 17.76, 18.20, 25.88). 
Since the Hamiltonian is dominated by the kinetic 
energy contribution Hq, all eigenvalues are close to these 
points shown in Fig. ^ 

B. ALS iteration methods 

We use two different quantities to monitor the con- 
vergence of the ALS iteration methods. As a necessary 
convergence condition for all iteration techniques the x n 
should tend to zero for n — ► oo. We therefore use the 
norm ||x„|| (maximal absolute column sum norm of x n ) 
as function of the number of iterations n as one conver- 
gence criterion. Thus, if this norm increases, the iteration 
method cannot converge. 

At each iteration step n the approximate solution 

n 

a n =^2xi (19) 

i=0 

is plugged into the decoupling equation (cf. Eq. ©) 

d n = a n PHQa n + a n PHP + QHQo n + QHP , (20) 

defining in this way the quantity d„. Its norm \\d n \\ rep- 
resents a deviation (distance) from the exact solution, 
where the decoupling equation is exactly fulfilled. This 
quantity should also tend to zero for n — > oo. As a second 
convergence criterion we calculate the eigenvalues at each 
iteration step by diagonalization which is surely possible 
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FIG. 2: The norm of x„ and of the deviation from zero d n 
of the decoupling equation for the model as function of the 
iteration number n using the standard ALS iteration method 
with one momentum cutoff. 



for the model and monitor their behavior as function of 
the number of iterations. 

We begin our analysis with the standard ALS iteration 
method as described in Sec. Ill Al In Fig. [21 the norm 
||a; n || and the deviation \\d n \\ from the exact solution of 
the decoupling equation are shown as a function of the 
iteration number n. After the first few iterations, these 
quantities increase with n signaling a divergence of the 
method but then decrease again for a larger number of 
iterations. This could lead to the incorrect conclusion 
that the ALS iteration method converges anyway to the 
right eigenvalues. 

This is demonstrated in the following: The eigenvalues 
in the P-space are given by the eigenvalues of p{a n ) and 
are plotted in Fig.|21as a function of the iteration number. 
The dashed lines are the final exact eigenvalues obtained 
by diagonalization of the model Hamiltonian. 

The starting values for the eigenvalues are given by 
PHP = p(Q). Due to our grid choice for the model we 
can identify the P- and Q-space eigenvalues uniquely. 
We violate in this way the convergence condition (II 1|) 
because the Q-space eigenvalue near the cutoff is smaller 
than the corresponding P-space eigenvalue (cf. Fig. 
Thus, we expect that the standard ALS iteration method 
does not converge. 

This is clearly visible in Fig. [3] for the eigenvalue Ai 
where the evolution of the P-space eigenvalue (filled 
squared) and the corresponding Q-space eigenvalue (open 
squared) are shown. Although the x„'s tend to zero, the 
eigenvalues converge for large n to the ones with min- 
imum absolute value. The P-space eigenvalue evolves 
from Ai ~ 17.78 to 15.05 and the Q-space one from 15.00 
to 17.76. This mixing between a P- and Q-space eigen- 
value during the decimation is typical for a divergence 
in the ALS iteration method. Due to our grid choice we 
have only one eigenvalue mixing the other three P-space 
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FIG. 3: The P-space eigenvalues evolution for the model using 
the standard ALS iteration with one momentum cutoff. The 
open squares label the nearest Q-space eigenvalue. See text 
for details. 



eigenvalues converge to the correct solutions. But in gen- 
eral an arbitrarily larger number of eigenvalues which vi- 
olate condition lllj l could mix. 

Unfortunately, for realistic interactions like the 
hyperon-nucleon interaction it is not possible to identify 
and sort the proper P- and Q-space eigenvalues for each 
ith-subblock. One has to find alternative criteria which 
signal an eigenvalue mixing. One alternative is the norms 
||.t„|| and ||d n ||. During the mixing the ||x n || and ||<i ra [| 
rise and once the eigenvalue sorting is completed they 
begin to drop again. Thus, a rising ||x n || (or ||e? n ||) and 
a following dropping with the iteration number n signals 
a level crossing and the standard ALS iteration method 
converges to wrong results. This behavior is shown in 
Fig.H 

As already mentioned, the standard ALS iteration 
method can be modified and improved by adding a diago- 
nal matrix M = diag(mi, . . .) to the decoupling equation 
with arbitrary chosen m; (cf. Eq. l|12|) . Then, the eigen- 
values converge closest to the chosen numbers mj. But 
the crucial point here is how to chose this matrix Ml No 
general criterion is known. For the model calculation it 
is possible to try several different ansatze, which indeed 
prevent the mixing of the P- and Q-space eigenvalues 
and result in stable eigenvalues. 

Unfortunately, all these ansatze cannot be applied to 
realistic interactions. For a realistic coupled channel cal- 
culation much more mesh points have to be taken into 
account. Then the identification of the P- and Q-space 
eigenvalues for each flavor channel i is not possible any 
longer. This is also visible in Fig. ^ O ne cannot dis- 
tinguish between the channel 2 P-space eigenvalues for 
small momenta k and the channel 1 P-space eigenvalues 
for momenta around the cutoff Ap. Even an alternative 
identification of the appropriate eigenvalues by means of 
eigenvector distributions, as described e.g. in Ref. |l4j . 
fails for realistic multi-flavor channel interactions. 
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FIG. 4: The kinetic energy Ho as a function of the momentum 
k for the model similar to Fig.Q The line with the squares de- 
notes the first-flavor channel, the other one the second-flavor 
channel. 



For our model the situation is different. By construc- 
tion, the eigenvalues of the P- and Q-spacc are well- 
separated and can be assigned uniquely. One way out of 
this dilemma is the introduction of a so-called channel- 
dependent cutoff. The use of the channel-dependent cut- 
off will cure not only the identification problem of the 
eigenvalues but also has a natural physical interpretation 
which we will discuss in the following. 



C. Channel-dependent cutoff 

We assume that the dominant contribution to the full 
Hamiltonian comes from its kinetic energy Hq and the in- 
teraction is a small perturbation. The channel-dependent 
cutoff is a general prescription how to define the model 
space for low-momentum interactions that include flavor 
transitions with different masses. This is for example of 
relevance for the hyperon-nucleon interaction in the par- 
ticle basis 0] . The idea of the channel-dependent cutoff 
is to truncate the kinetic energy of each channel to the 
same energy Ap. The truncation of the kinetic energy 
in this way leads to different momentum cutoffs in each 
flavor channel Aj . This is depicted in Fig. 0] 

In principle there are two possibilities how to construct 
the momentum cutoffs. One possibility is to fix the value 
of the kinetic energy Ap and then to calculate the corre- 
sponding momentum cutoffs Aj. But the following pro- 
cedure is better adapted to a generalization of the V\ OV! k 
method to multi-flavor channels. Instead of fixing the ki- 
netic energy a priori, one calculates its value for a given 
momentum cutoff of the flavor channel with the smallest 
mass. For our model this means to choose Ai for channel 
1 and calculate A2 for channel 2 by requiring 

#0,1 (Ai) - ffo, 2 (A 2 ) . (21) 
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For example, choosing a momentum cutoff Ai = A = 5 
for channel 1 this results in a smaller cutoff in channel 2 



A 2 = ^2^ 2 ^—L+ mi - m2 j « 3.1, (22) 

which is denoted by the vertical line in Fig. 21 Due to 
the new momentum cutoff A 2 the mesh points have also 
changed for channel 2. A generalization to multi-flavor 
channels is straightforward. 

Using this channel-dependent cutoff, the eigenvalues 
of the P- and Q _s Pace are now clearly separated and the 
necessary convergence conditions and l|14fl are ful- 
filled. Here, we have assumed that the potential is a 
perturbation and does not influence the eigenvalue sort- 
ing strongly. We have checked this assumption by a nu- 
merical calculation of the eigenvalues also in the realistic 
calculation. 

As already mentioned, note that with only one momen- 
tum cutoff for multi-flavor channels it is not possible to 
fulfill these convergence conditions and the ALS iteration 
methods must therefore fail. 

Another, more physical motivation for the channel- 
dependent cutoff is the following: The derivation of 
the Mow k renormalization group equation is based on 
a Lippmann-Schwingcr equation for the T-matrix. The 
propagator in the LSE is inverse proportional to the en- 
ergy, which is truncated by a momentum cutoff. For 
channels with different reduced masses this would allow 
for different kinetic energies. Thus, it is natural to re- 
strict the kinetic energy available in the propagator and 
not to restrict the momenta. For one-flavor channels this 
difference does not play any role. 

This becomes more transparent if one considers an on- 
shell scattering process of particles with different masses. 
Assume the scattering of a particle of channel 2 into a 
particle of channel 1 with a relative incoming momentum 
around the model space cutoff, i.e. k' w A and masses 
m 2 > mi. Since we consider on-shell scattering, the en- 
ergy 

A 2 

ff , 2 (fc'«A)« — +m 2 = E. (23) 

2/U2 

is fixed, say to E. This would result in 

E = JTo,i(*) = T- + mi. ( 24 ) 

for the scattered particle 1, which in turn yields for the 
outgoing relative momentum 



= y y«i (j- + 2(m 2 - mi)j . (25) 

Since we have assumed m 2 > mi this outgoing rela- 
tive momentum is larger than our cutoff A. For example, 
choosing for the incoming momentum k' = 4.9 the out- 
going momentum would be k ~ 7.4. Since all momenta 
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FIG. 5: Same as Fig.[5]for the model using the standard ALS 
iteration method with a channel-dependent cutoff. 



are limited by the cutoff A = 5 such transitions are not 
allowed. This is also clearly visible in Fig. 0] 

Since the channel-dependent cutoff restricts the kinetic 
energy it can be represented by a horizontal line in this 
figure. One never crosses this line in any on-shell scat- 
tering process. 

When we apply the channel-dependent cutoff to 
the model, the mesh points will change due to the 
different momentum cutoff, which is shown in Fig. 0] 
This yields also a change of the eight exact eigen- 
values. In this case we find the eigenvalues H = 
diag(9.55, 14.50, 15.05, 24.04, 12.96, 14.56, 14.88, 20.22) 
which arc again close to the points in Fig. 

Due to the channel-dependent cutoff the convergence 
criterion for the ALS iteration methods are fulfilled now 
because the P-space eigenvalues are these eigenvalues 
which arc closest to zero. The different iteration schemes 
converge monotonic and to the correct P-spacc eigenval- 
ues. An example for the standard ALS iteration is shown 
in Fig. \5\ 

The evolution of the four P-space eigenvalues for the 
model as a function of the number of iterations is dis- 
played in Fig. EJ The dashed lines are the exact eigen- 
values. All eigenvalues converge to the correct ones. Due 
to the small gap between the largest P-space and the 
smallest Q-space eigenvalue in the model ,the number of 
iterations of is slightly larger compared e.g. to Fig. 



IV. SUMMARY 

In this work the ALS iteration methods and its mod- 
ifications are investigated, which are necessary to solve 
coupled-channel equations with different masses. Instead 
of working with realistic interaction potentials, a simple 
and solvable model, which mimics all the relevant prop- 
erties of a realistic interaction is introduced. 

First, we applied the ALS iteration method and its 
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FIG. 6: The P-space eigenvalue evolution for the model sim- 
ilar to Fig. 01 using the standard ALS iteration method with 
a channel-dependent cutoff. 



modifications to the one-flavor channel and discussed its 
failure for multi-flavor channels. Due to the different 
masses inherent in the multi-flavor case the necessary 
convergence conditions are violated. 

By an introduction of a channel-dependent cutoff we 
could generalize these ALS iteration modifications and 
could restore the convergence. This channel-dependent 



cutoff results in an effective interaction which is still 
energy-independent and thus preserves a major advan- 
tage of the Mow k potential. 

An on-shell scattering example of two particle with 
different masses provides a more physical motivation for 
this novel cutoff. Without the channel-dependent cutoff 
unphysical scattering processes out of the model space 
become possible. First encouraging results for a realis- 
tic hyperon-nucleon interaction by means of the channel- 
dependent cutoff have been published in Ref. Fur- 
ther calculations including a broader variety of input po- 
tentials such as the Nijmegen potentials (NSC97 jllj . 
NSC89 0), Jiilich potential (Juelich04 and po- 

tentials based on chiral symmetry (LO potential 
will be presented elsewhere |Tfjj . 
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